Thermodynamics of Solitonic Matter Waves in a Toroidal Trap 



O 
O 

(N 



5h 



C3 



S3 
O 

o 



> 

o 

\o 
o 

-I— > 



i 

o 
o 



X 
S3 



L. Salasnich 1 , A. Parola 2 , and L. Reatto 3 
1 CNR-INFM and CNISM, Unita di Milano Universitd, Via Celoria 16, 20133 Milano, Italy 
2 Dipartimento di Fisica e Matematica, Universitd dell'Insubria, 
Via Valleggio 11, 22100 Como, Italy 
3 Dipartimento di Fisica and CNISM, Universitd di Milano, 
Via Celoria 16, 20133 Milano, Italy 

We investigate the thermodynamic properties of a Bose-Einstein condensate with negative scat- 
tering length confined in a toroidal trapping potential. By numerically solving the coupled Gross- 
Pitaevskii and Bogoliubov-de Gennes equations, we study the phase transition from the uniform 
state to the symmetry-breaking state characterized by a bright-soliton condensate and a localized 
thermal cloud. In the localized regime three states with a finite condensate fraction are present: 
the thermodynamically stable localized state, a metastable localized state and also a metastable 
uniform state. Remarkably, the presence of the stable localized state strongly increases the critical 
temperature of Bose-Einstein condensation. 



PACS Numbers: 03.75.Kk 

In recent experiments a repulsive Bose-Einstein con- 
densate (BEC) has been produced and studied in a quasi 
one-dimensional (ID) ring [1,2]. These experiments and 
also previous experimental investigations with cold atoms 
in torodial traps [3] are important steps to achieve ultra- 
high precision sensors with atom interferometry (see for 
instance [4]). The case of an attractive BEC in a ring 
has not yet been experimentally investigated but ap- 
pears very interesting: a quantum phase transition from 
an azimuthally uniform condensate to a bright soliton 
has been predicted [5,6]. Dynamically stable multi-peak 
bright solitons may also appear [7]. 

In this paper we investigate the effect of tempera- 
ture on the uniform-to-localized phase transition by solv- 
ing the Gross-Pitaevskii equation (GPE) of the macro- 
scopic wave function of the Bose condensate and the 
Bogoliubov-de Gennes equations for the quantum deple- 
tion and the thermal cloud. We show that, for a fixed 
number TV of atoms, the uniform solution always exists 
but it is thermodynamically stable only up to a critical 
interaction strength which depends on the temperature. 
Above this critical interaction strength both condensate 
and thermal cloud become localized. In addition, we de- 
termine the finite-temperature phase-diagram of the sys- 
tem for a sample of alkali-metal atoms in a ring. Our 
results are relevant not only for the physics of cold atoms 
but also for the nonlinear science of solitons. In fact, a 
mixture of condensed and non condensed attractive par- 
ticles [8] correspond to incoherent solitons of light in Kerr 
media [9]. 

We consider a Bose gas with negative scattering length 
(a s < 0) confined in a toroidal potential and model 
the transverse confinement with a harmonic potential 
of frequency oj±. The two characteristic lengths of the 
toroidal trap are the azimuthal radius R and the trans- 
verse harmonic length a±_ = {%/ 1 (muj^)) 1 / 2 . To avoid the 
confinement-induced resonance at \a s \ ~ ai [10] we im- 
pose that \a s \ <C a± [11]. At low temperature T and with 



a finite number N of atoms, in the Bose gas there is the 
thermal cloud but also the BEC. Under the condition 
R ^> a± the azimuthal wavefunction <f>(z) of the BEC 
in the ring satisfies the ID non-polynomial Schrodinger 
equation (NPSE) [7,12] 
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where z — RB is the azimuthal coordinate with 9 the 
azimuthal angle. This equation has been deduced from 
the 3D GPE by using a Gaussian transverse wave func- 
tion with a density-dependent width and neglecting cur- 
vature effects [7,12]. The non-polynomial term fj,(no(z)) 
is a function of the BEC density tiq(z) — A^ 1 ^ ( ^ ) 1 2 and it 
is given by fi = dS/dno, where £ = hcj± no(l + 2a s no) 1 ^ 2 
and a s < is the 3D s-wave scattering length [7]. The 
wave function <p(z) is normalized to one and satisfies the 
condition of periodicity 4>(z + L) = 4>{z), where L = 2ttR. 
The chemical potential /2 is fixed by the normalization 
condition J \<fi(z)\ 2 dz = 1. Quantum and thermal deple- 
tions of the BEC are obtained by solving the Bogoliubov- 
de Gennes (BdG) equations [13] for the quasi-particle 
amplitudes u a (z) and v a (z), given by 
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where e a are the energies of quasi-particle excitations and 
the operator C is 
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The local density of the quantum depletion reads 
n m i.t{z) = Y] a \v a (z)\ 2 , while the thermal local density 
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is n T (z) = J2a {\ u <x{z)\ 2 + \Va(z) 

(exp(e Q / ksT) — 1) 1 are the Bose numbers of occupa- 
tion for the non-interacting quasi-particles, ks is the 
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Boltzmann constant and T is the absolute temperature 
of the thermal cloud. The number N of non-condensed 
atoms is thus given by N = N out + Nt, where N out — 
J n out {z)dz and Nt — J nr{z)dz. In the weak-coupling 
limit, where a s no <C 1, the energy density E becomes 
£ = hui±n + hui±a s nl, the NPSE (1) reduces to the fa- 
miliar ID GPE [7] and the BdG equations (2,3) reduce 
to the BdG equations of the strictly ID problem. Note 
that the Eqs. (1) and (2,3) could be improved by taking 
into account Popov and Bcliaev corrections, which are 
usually not negligible only when the condensed fraction 
is very small [14]. 

In our approach, by fixing the total number N — 
No + N of atoms and the temperature T, the equilib- 
rium configurations are found from the Hclmholtz free 
energy [13] 



F = E + E + p(N - N ) - TS 



where 
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= N J 



(4) 



dz (5) 



is the BEC energy. The out-of-condensate energy E and 
entropy S are instead given by 



n a - / \v a (z)\ 2 dz 
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S=-k B J2 [n a log(n a ) - (1 + n a ) log(l + n a )\ . (7) 

The solution of the ID GPE with periodic boundary 
conditions is known analytically [5]: it is uniform for 
< 7./V0 < ir 2 a±/L, with 7 = 2|a s |/a^ the interatomic 
strength, and becomes localized for 7A0 > ir 2 a±/L. By 
using the NPSE we find that the transition strength ap- 
proaches the ID GPE one for large L [7]. The main 
difference between ID GPE and NPSE is that the NPSE 
correctly gives the collapse of the localized solution for 
jN > 4/3, while the ID GPE docs not predict col- 
lapse [7]. It is important to observe that the uniform- 
to-localized transition depends on iVo and only implicitly 
on N — No + N out + Nt- We shall show that, for a fixed 
number ./V of atoms, the uniform solution exists for any 
value of the interatomic strength 7 but it is thermody- 
namically stable only below a critical strength j c . Above 
7 C the stable state is a localized solution that minimizes 
the free energy. 

The BdG equations (2,3) are easily solved if </>(z) is 
uniform: 4>(z) = 1/vT. In this case the quasi-particle 



amplitudes Uk{z) — u^e lkz and Vk(z) 
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are plane 



waves and the wave vector k is quantized: k — (2tt/ L)j, 
with j an integer number. When <f>(z) is not uniform, 
one must numerically solve Eq. (1) and Eqs. (2,3). We 
use a finite-difference space discretization and diagonalizc 
the matrix associated to the operator C of Eq. (3). In 



our calculations we use various matrix dimensions up to 
4000 x 4000. The numerical procedure is discussed in 
[15]. In our approach the transverse excitations do not 
contribute to the thermodynamics: this assumption is 
reasonable only if the system is quasi-lD. 
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FIG. 1 Free energy F as a function of the scaled inter- 
atomic strength 7 = 2\a s \/a± for an attractive Bose gas 
of N = 10 4 atoms in a toroidal trap with L/a±_ = 100. 
Dashed line: uniform solution. Solid line: localized solu- 
tions. The energies F and ksT are in units of Huj±. 

In Fig. 1 we plot the results of our numerical calcula- 
tions of the free energy F as a function of the scaled inter- 
atomic strength 7 = 2|a s |/aj_ for increasing values of the 
temperature T. We choose N = 10 4 atoms and a toroidal 
geometry with L/a± — 100. While the uniform solution 
(dashed line) exists for any strength 7, the localized so- 
lution (solid line) exists only above a critical strength j e . 
As expected, at T ~ we find 7 C ~ n 2 a±/L = 0.0987. 
For a given 7 (7 > 7 e ) there are two localized states: 
one is thermodynamically stable (lower solid line) and 
the other is metastable (upper solid line). The metasta- 
bility is guaranteed by the Bogoliubov eigenvalues being 
real. At very-low temperature, see panel (a) of Fig. 1, 
for 7 > 7 e the solid line of the metastable localized so- 
lution is very close to the dashed line of the metastable 
uniform solution. At higher temperature, see panels (b) 
and (c) of Fig. 1, the free energy of the metastable local- 
ized solution, but also the free energy of the stable one 
at small 7 > "f e , is higher than the free energy of the uni- 
form solution. The strength -f c at which the dashed line 
crosses the lower solid line determines the point where the 
first-order transition takes place. At T ~ we find that 
7 C ~ 7 e . The strength 7 C grows with the temperature T 
but in general 7 C does not coincide with the strength 7 e at 
which the localized solutions appear. Thus, at finite tem- 
perature, for < 7 < 7c there is only the stable uniform 
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state; for 7 e < 7 < 7c there are two metastable localized 
states and the stable uniform state; for 7 > -f c the sta- 
ble state is one of the two localized state while the other 
localized state and the uniform state arc metastable. 
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FIG. 2 Condensate fraction No/N as a function of the 
scaled interatomic strength 7 = 2\a s \/a± for N = 10 4 
atoms with L/a± — 100. Dashed line: uniform solution. 
Solid line: localized solutions. The vertical dot-dashed 
line indicates the strength N^ c at which the uniform-to- 
localized phase transition takes place. 

A coherent bright soliton and a thermal cloud coexist 
in the stable state, as in the metastable localized state. 
This is clearly illustrated in Fig. 2 where we plot the 
condensate fraction Nq /N as a function of 7 for different 
temperatures. At low temperature, see panel (a) of Fig. 
2, the uniform state (dashed line) is a quasi-pure BEC 
for < 7 < 7 e . For 7 > 7 e the condensate fraction of the 
uniform state quickly decreases. At low temperature the 
uniform-to-localized phase transition does not produce 
a relevant quantum depletion N out /N at the transition 
point. This result is obtained with a fixed number of 
atoms (N = 10 4 ). We have verified that by fixing Nq in- 
stead of N (as done in [5]) the condensate fraction N /N 
goes to zero for 7A0 ~ n 2 a±/L. 

At higher temperature, see panels (b) and (c) of Fig. 
2, BEC and non-condensate component coexist also in 
the uniform state. Fig. 2 shows that the condensate 
fraction of the localized states is higher than the con- 
densate fraction of the uniform state. Remarkably, at 
the transition strength j c there is a sizeable jump in the 
condensate fraction. Moreover, for the stable localized 
state the ratio N /N initially grows by increasing the 
interatomic strength. As previously discussed, the BEC 
soliton of the localized state will collapse at 7A0 — 4/3, 



but the true collapse is avoided by populating the out-of- 
condensate cloud with N = N — N atoms for A 7 > 4/3. 
We have verified that in the stable localized state the 
number No of condensed atoms is always very close to 
the collapse value 4/(37) an d, for a fixed 7, the number 
N has a very small temperature dependence. It follows 
that one can estimate the condensed fraction of the sta- 
ble localized state as N /N ~ 4/(3A7). This estimation 
is not so rough: choosing for instance A7 = 2 one finds 
No/N ~ 2/3 in good agreement with the result shown in 
the upper solid line of Fig. 2 (see where N /N — 0.66). 

In the previous calculations, transverse modes have 
been neglected and the thermally excited quasi-particles 
are allowed to populate only the lowest transverse state. 
This assumption can be justified by a simple argument: 
In the weakly localized regime (low Nj) the spectrum 
of elementary excitations is well approximated by the 
free Bose gas, where the logitudinal modes have energies 
0(hu!±a 2 L / L 2 ) much lower than the O(hoj±) transverse 
modes. In this limit we have explicitly verified that, with 
the chosen parameters (L/a± = 100), for fc^T < 5hui± 
the number of quasiparticles in the excited transverse 
modes is indeed negligible. In the strongly localized re- 
gion, both longitudinal and transverse low energy excita- 
tions are known in the large L limit [6,17], where a finite 
0(huj±) gap separates the two branches. Moreover, our 
calculations show that for A 7 < 4/3, thermal depletion 
does not affect deeply the thermodynamics of the model 
up to the transition temperature. 
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FIG. 3. Azimuthal density profiles of the thcrmody- 
namically stable state for N = 10 4 atoms with L/a± = 
100 and kBT/(Tiuj±) = 2. Length z is in units of 

a± = (VK)) 1/2 . 

For completeness, we have repeated the calculations by 
using the ID GPE. Also in this case we have found two 
localized solutions but the stable one always maintains 
a condensate fraction close to one, due to the fact that 
the ID GPE does not predict the collapse. Thus the 
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transverse structure of the soliton, taken into account by 
the NPSE, strongly modifies the thermodynamics of the 
attractive Bose gas. 

In Fig. 3 we plot the density profile no(z) of the con- 
densate cloud (solid line) and of the density profile tit(z) 
of the thermal cloud (dot-dashed line) for the stable local- 
ized state and two values of 7. We have verified that the 
density profile n out {z) of the quantum depletion reduces 
by increasing the temperature T and that the quantum 
depletion N out /N becomes relevant only when the ksT 
is of the order of the lowest Bogoliubov energy level. 
This energy can be easily estimated as ksT/ (hw±) = 
(h 2 /{2m))(2ir/L) 2 /(huj ± ) = 2ir 2 a\/L 2 ~ 0.002. 



that our predictions can be verified by using a mag- 
netic ring with geometric parameters not too far from 
those used in Ref. [1] for trapping a repulsive conden- 
sate of 87 Rb atoms. For instance, a sample of 7 Li atoms, 
which have a negative scattering length (a s = — 1.4 nm), 
can be trapped in a ring with transverse width a± ~ 6 
/im, which gives huj^/ks = Ti 2 /(2a\) ~ 2 nK. In this 
way, having TV = 10 4 atoms the interaction strength is 
A^7 = N2\a s \/a± ~ 5. Then, by taking the azimuthal 
radius R = L/(2tt) ~ 100 nm one can test the predic- 
tions of Fig. 4 varying the temperature or tuning the 
scattering length a s around a Feshbach resonance. 
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FIG. 4 (color online). Phase diagram in the plane tem- 
perature T vs interatomic strength 7 for N = 10 4 atoms 
in a toroidal trap, with L = 100. The energy ksT is in 
units of huj± and the length in units of a± . 

In Fig. 4 we plot the phase diagram of the atomic 
cloud of bosons in the plane (JV7, T) for L/a± = 100 
and N = 10 4 . In the figure we insert a dashed line 
at ksT/ (hiu±) — 6.01, which is the transition temper- 
ature Tbec of the Bose-Einstein condensation above 
which the condensate fraction of the uniform state is zero. 
The transition temperature Tbec reduces by increasing 
L/a±. In fact, it is easy to show [16] that ksTBEc/(^±) 
scales as Na]_/L 2 . Fig. 4 displays as a solid line the curve 
of the uniform-to-localized transition. Remarkably, the 
BEC transition temperature of the localized solution is 
much larger than the uniform one. Fig. 4 shows that 
with A^7 > 0.68, there is a first order transition from 
localized BEC phase to the the uniform thermal phase 
without BEC. Instead, for N7 < 0.68 the no-BEC phase 
is obtained starting from the uniform phase by increasing 
the temperature. 

In conclusion, we have shown that an attractive Bose 
condensate in a quasi one-dimensional toroidal trap dis- 
plays at finite temperature new and intersting features, 
like the coexistence of a coherent self-bound bright soli- 
ton with a thermal cloud, the avoiding of true collapse 
via population of the out-of-condensate (thermal) com- 
ponent, and the enhancement of the Bose-Einstein tran- 
sition temperature for the localized solution. We stress 
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